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1. Introduction 



The c = 1 string (an excellent review is Ref. |T[]) is a perturbatively consistent string 
theory in two spacetime dimensions. One of its attractive features is that it is 
solvable: from the powerful techniques of Matrix Quantum Mechanics (MQM), cor- 
relation functions of the momentum modes ("tachyons") can be determined to all 
orders in the string coupling (inverse cosmological constant). This holds true even 
in the Euclidean theory at finite radius R. Another feature is that its generalisation 
to the type noncritical string, has similar properties in perturbation theory but is 
believed to also be non-perturbatively well-defined. 

This makes the c = 1 string and its cousins a good laboratory to study var- 
ious open questions in string theory. Two such questions that we would like to 
understand better in the noncritical context are the properties of string-scale black 
holes, and the nature of various dualities, including open-closed string duality || [|]. 
Much work has been done on the former (some interesting recent studies can be 
found in Refs. || [[/J ) , while the latter question has also yielded some important 
illuminations @ |2j @ [|| . 

It is known |TjJ[jnj that basic properties of black holes in noncritical string the- 
ory are controlled by condensates of winding tachyons in the Euclidean-continued 
background. These are thermal tachyons: strings winding around the compact time 
direction. It would therefore be useful to know the correlators of winding modes 
in Euclidean noncritical string theory to all orders in the string coupling (and even 
nonperturbatively in the stable type-0 case) as a function of /i and R, where /i -1 is 
the inverse string coupling and R is the radius of the Euclidean direction (inverse 
temperature). From the matrix model point of view, winding modes are related to 
the nonsinglet sector of the model, in which the eigenvalue fermions are no longer free 
but mutually coupled [fL7||||18| . Computing correlators in this way is a harder task[19 



and has raised some new puzzles involving leg factors which we will discuss in a later 
section. But one way to find the desired correlators is to assume that T-duality holds 
and perform it on the momentum correlators. This provides one of our motivations 
to study momentum correlators in the Euclidean theory in more explicit detail than 
has already been done. 

As mentioned above, momentum correlators in the Euclidean c = 1 string are 
known in principle. They are summarised in the Toda hierarchy or symmetries |2"U 
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or Hirota bilinear equations, or Normal Matrix Model (NMM)|2T[, all of which are 
supposed to be mutually equivalent. For the special case of self-dual radius R = 1 
of the Euclidean time direction, they are encoded in a Kontsevich-Penner matrix 
model|22|] [^J (see also P^||[25|| ). We will summarise some relevant information about 
these solutions below. But while all these formal solutions allow us to extract the 
perturbation series for any specific correlator after a sufficient amount of work, we 
do not have many explicit answers in terms of special functions depending on the 
radius R and inverse string coupling \i. 

At finite radius, correlators have been computed mostly at tree-level (corre- 
sponding to the dispersionless limit of the Toda hierarchy) or to a few low orders 
in perturbation theory. For example, while the 2n-point function of n unit winding 
modes and n anti-winding modes is known as a function of n and R at tree level 
[f2(J[|n|, an explicit expression for the same correlator to all orders in perturbation 
theory does not seem to exist in the literature 1 . To be more specific, denote by T q 
the tachyons of momentum q = n/R, and by T q the tachyons of n units of winding, 
where q = nR is the value ofpi = —Pr in vertex operator language. An explicit form 
is known for ((T_i/^)"(Ti//j) n ) at tree level. The T-dual of this expression was used 
in Ref. [0 to extract the critical behaviour of the Sine-Liouville theory defined by 
perturbing the original c = 1 string with T_r + 7r and then tuning the cosmological 



constant fi to zero. In particular, Ref.|16| showed that a sensible theory exists after 
this tuning, but only when the radius of the Euclidean direction lies in the range 
1 < R< 2. 

One would like to know the structure of this correlator to all string loop orders. 
Accordingly, in what follows we will study ((T^i/R) n (Ti/R) n ) in detail, and one of 
our main results will be a simple formula for this correlator as a function of fi and R 
for every n. We expect this to lead to a better understanding of the exponentiated 
correlator (exp(T_x/ R + Tya)), which in turn is T-dual to the vortex condensate 
(exp(71/j + Tr)) that relates directly to Euclidean 2d black holes. 

Another motivation for our work is to understand T-duality of the c = 1 matrix 
quantum mechanics. This is established at the level of spectrum of states, since 



the partition function without perturbations is known to be T-dual |29]. Also, a 



X A differential equation for these correlators was written down injl6| together with an iterative 
solution to a few orders. Related work on Euclidean correlators can be found in Refs.|27|| JE8|. 
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formal argument has been given [16] that the winding correlators, like the momentum 
correlators, are given by a Toda hierarchy 2 . However, to our knowledge, beyond this 
result and a computation in |19| , there has been no direct comparison of correlators in 
the momentum and winding sectors 3 . A convincing test of T-duality would consist 
of computing pure-momentum correlators in terms of free fermion eigenvalues, T- 
dualising the answers and comparing them with pure-winding correlators computed 
from the nonsinglet Hamiltonian. Ideally this should even be done beyond tree 
level. Although we will not be able to carry out such a test here, we have tried to 
systematise one side of the duality in a way that can be eventually compared with 
the other side when nonsinglet computations become more practicable. 

In particular, the most direct way to check T-duality comes from comparing 
two-point functions. Accordingly we work out all two-point functions of momentum 
modes. In Ref.|31J, the two-point function of unit-momentum modes was computed 
and an attempt made to match the leading result with a computation in the first 
nonsinglet sector of the matrix model, namely the adjoint sector. The compari- 
son revealed the presence of unexplained normalisation factors. It was pointed out 
in Ref.0] that if one could compute two-point functions of more general winding 
modes, namely (T-nR%iR), one might be able to shed some light on these normalisa- 
tion factors. With this motivation we have performed this computation and obtained 
a simple explicit result, again as a function of \i and R and for all n. In a later section 
we discuss the relation to the non-singlet sectors. 

Our initial computations have been performed using both the MQM and a model 
of constant matrices called the Normal Matrix Model (NMM) []2T| , with perfect agree- 
ment between the answers. In the former case we used the known infinite-radius 
correlators in the physical MQM (real and noncompact time) |JT| , and a formula 



which converts these to the correlators for Euclidean compact time||32||. In the lat- 
ter case, we will describe how one performs computations after reducing the NMM 
to eigenvalues. One surprise emerging from comparison of the two approaches is 
that the NMM successfully computes correlators even when the matrices are of finite 



2 For a discussion of T-duality in type 0A,B matrix models, see Ref . ]30[ . 

3 In addition to pure momentum or pure winding correlators, one would also like to know the 
correlators for a mixture of momentum and winding modes. In this case one has no choice but to 
tackle the difficult nonsinglet sector problem. The system is not expected to be integrable and the 
correlation functions are not known so far. 
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rank N, a stronger property than was claimed in Ref.[21], who did however suggest 
that the model contains some information even at finite N. We find that it actually 
contains complete information at finite N in the following sense: given a correla- 
tor, there is a minimum value N min such that this correlator when computed in the 
NMM gives the correct result, to all orders in l//i 2 , for all N > N m i n . This makes 
the NMM a potentially powerful combinatoric tool. We then go on to demonstrate 
its power by deriving a combinatoric formula for the general correlation function 
((T. 1/R ) n (T 1/R ) n ) for any n. 

We start in Section |2] by describing the two relevant matrix models, Matrix 
Quantum Mechanics and Normal Matrix Model. The former is too well-known to 
need a detailed discussion and we skip directly to the calculational techniques and 
answers. For the latter, we review the model in some detail, with special attention to 
the role of the matrix rank N. In Section || we work out some relevant correlators as a 
function of \x and R from MQM. In Section |] we reproduce these correlators from the 
NMM, where we note the phenomenon that for a fixed correlator, the NMM at any iV 
greater than a minimum value gives the complete answer. After a discussion of why 
this works, we use this property to derive a combinatoric formula for correlators of 
any number of unit momentum modes. In Section |5| we discuss applications of these 
results to some physically interesting problems, and conclude in Section || Several 
calculational details are presented in the appendices. 

2. Matrix Quantum Mechanics and Normal Matrix Model 
2.1 Matrix Quantum Mechanics 

Matrix Quantum Mechanics is a model of a single N x N hermitian time-dependent 
matrix M(t). In the absence of perturbations, the partition function of the model is 
given by: 



Zmqm = J[dM}exp -N J dt tr ((D t M) 2 + M 2 ) 



(2.1) 



where D t M = M + zL4.£,M] is the covariant derivative with respect to the time 
component of a gauge field. 

The gauge field acts as a Lagrange multiplier and projects the model to the singlet 
sector, which is a system of N non-interacting non-relativistic fermions moving in 
an inverted harmonic oscillator potential. In the double-scaling limit, the fermi sea 



- 5 - 



is filled nearly to the top and the number of fermions is taken to infinity. The 
scaled distance to the top of the potential, /i, is kept finite and corresponds to the 
cosmo logical constant. This model provides a description of 2D string theory, with 
playing the role of the string coupling g s . 

The physical modes of 2D string theory can be constructed in terms of fermion 
eigenvalues. In |H] this model was used to calculate correlation functions of c = 1 
string theory at infinite radius. One starts by computing correlators of free-fermion 
bilinears, which in turn can be used to extract correlators of the loop operators: 

0(k,£) = J dte lkt tre- m{t) (2.2) 

Extracting the leading behaviour of these loops for small £, one has 

0(k,£) ~£ |fc| T fc (2.3) 

The Tfc are identified with the c = 1 string theory tachyons. When compared with 
the corresponding operators in Liouville theory, there is a change of normalisation: 

Tk\MQM = T(l^l) Tk | Liouville (2.4) 

However this fact will not be relevant for us, since in what follows we will always 
work with the operators Tk in the MQM basis, i.e. the LHS of the above equation. 
When the time direction is Euclidean and compact, we are in the finite tem- 



perature theory. Starting from the infinite-radius correlator, one can show [32] that 
correlators in the Euclidean theory at finite radius are obtained as: 

(T qi T q2 ■ ■ ■ Tq n )n = — -. . (T qi T q2 ■ ■ ■ TqJ)^ (2-5) 
sin 

In addition one must replace the momentum-conserving 5-function as: 

*(£>) ^R5 Eiqh0 (2.6) 

i 

The above prescriptions follow from the fact that the compact radial direction 
introduces an additional factor in the loop momentum integrals of the infinite-radius 
calculation, and this factor can now be taken out of the integrals whence it becomes 
a differential operator acting on the infinite-radius answer 4 . 

4 It is also possible to calculate correlators directly at finite radius using the "reflection coefficient" 
formalism of Ref. |2(]] . Though we will not use this here, it would be interesting to know if our explicit 
results follow as easily in that approach. 
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In the finite-temperature theory, the above modes can be thought of as carrying 
"momentum" in the time direction. In this situation one also expects to find winding 
modes corresponding to the thermal scalars of finite-temperature string theory. Many 
physical properties of string theory are encoded in these degrees of freedom, which 
are therefore quite important to study. To find them in the matrix model we must 
go beyond the singlet sector, in which the gauge field is topologically trivial and can 
be gauged away. Consider the gauge-invariant Wilson-Polyakov loop variable: 



(2.7) 



where the trace is performed in the representation 1Z of SU(N). When 1Z is the 
fundamental representation, this is to be associated with a unit winding mode: 

W n =N ~ T R (2.8) 

Similarly the trace in the anti-fundamental will be 7L#. One can also have loops 
where the trace continues to be in the fundamental but the contour winds multiple 
times over the Euclidean time direction. Computation of the correlation functions 
of all these Wilson-Polyakov loops is done by observing that in their presence, the 
matrix model receives contributions from definite non-singlet sectors. In these sectors 
it reduces to eigenvalue fermions but now with mutual interactions. For example, 
the two-point function of unit winding modes can be identified as follows: 

(T- R T R ) ~(W R W N ) ~(W adjoint > (2.9) 



Liouville theory 



MQM 



MQM 



Thus computing the partition function of MQM in the adjoint sector determines 
the two-point function of winding modes. Since in principle this is an independent 
computation from that of the momentum tachyon correlators, it can actually be used 
to check T-duality of the c = 1 string. We will return to this issue in a subsequent 
section. 



2.2 Normal Matrix Model 

The Normal Matrix Model (NMM)|21[] is a relatively simple model of a complex 
matrix Z and its Hermitian adjoint, with the constraint that the two commute (hence 
Z is said to be "normal"). The potential is polynomial with an additional logarithmic 
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piece. The matrix Z is constant rather than time-dependent, so in this sense it is 
more similar to the c < 1 string backgrounds which do not have a time direction 5 . 

The NMM is proposed to describe the correlators of the c = 1 string to all orders 
in perturbation theory, as follows. Let us introduce its partition function: 

(2.10) 

Here R, v are some (in general, complex) parameters, which will correspond to the 
compactification radius of Euclidean time and the cosmological constant respectively. 
The parameters t k , t k are couplings to the gauge-invariant operators trZ fc , tiZ^ k and 
Z, Z^ are N x N matrices satisfying: 

[Z,Z*]=Q (2.11) 

The operators trZ k , trZ^ are identified with the tachyons T k / R , T_ k / R of momentum 
±-j| respectively. 

Since the matrix Z commutes with its adjoint, the two can be simultaneously 
diagonalised. The diagonalising matrices drop out of the action leaving behind Van- 
dermonde factors. It turns out that one gets a single power of the Vandermonde for 
the eigenvalues z±, Z2, ■ ■ ■ , zn of Z, together with its complex conjugate corresponding 
to Z^. Thus, for example, the partition function at t k = t k = is: 

N 

\{d 2 z t J] \z t - z 3 \ 2 e -^f =1 (^^)«+(^-iv + %i)Ef =1 iog^^ (2 . 12) 

i=l i<j 



■'NMM 



with an obvious generalisation to include the tachyon perturbations. 

At t k — t k — 0, it can be shown (though not directly from the action) that the 
NMM is invariant under the T-duality operation: 

R^^, fi^fiR (2.13) 
R 

This invariance is broken by the presence of momentum modes. Indeed, after T- 
duality, the tachyons T± k / R of the c = 1 string turn into winding modes of ±k units 
of winding, or equivalently (in vertex-operator language) of left /right momentum 
(Pl,Pr) — ±(kR, —kR). In what follows, these modes will be denoted T kR) T_ kR . 

5 Perhaps this is the underlying reason why the NMM describes Euclidean c = 1 strings at an 
arbitrary radius R, but does not have a simple R — > oo limit where one might recover the Lorentzian 
theory. 
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In |H] two distinct equivalences between the NMM and the c = 1 string were 
proposed. The first, referred to as "Model I", requires us to take the large- N limit 
of the NMM. The result in this case was that: 



Z c=1 GM,t) = lim 4mm(^M). (2-14) 

after the analytic continuation v = —ifi. 

However, another equivalence, "Model II", was proposed which did not involve 
a large- N limit. It was argued that the c = 1 string theory can be obtained from the 
NMM at finite N, provided v is set to the special value ^ (note that this corresponds 
to an imaginary cosmological constant): 

Z c= \ (/i = i-ft,t, t) = Z NMM {y = j[,t,t) , (2.15) 

In other words, the claim 6 is that an NMM calculation for a fixed integer value of N 
determines Z c= \ for a particular (imaginary) value of /i, namely 

N 

H = i- (2.16) 

If we T-dualise the above considerations so that t,t become couplings to winding 
tachyons, this relation becomes 

li = iN (2.17) 

The above results seem to indicate that for finite N we can only generate the answer 
at a fixed fi, in which case we would never obtain the perturbative expansion in 
powers of l//x 2 . However, below we will compute winding correlators using the 
NMM, and will see that it turns out much more powerful than expected. It actually 
does reproduce the entire perturbative correlators, as functions of \i and R, even at 
finite values of N. Evidence for this fact, as well as an explanation of it, will be 
provided in subsequent sections. 



6 The authors of Rcf. stated this a little differently: that one obtains c — I string amplitudes 
as a function of fj, by computing NMM correlators as a function of N and fi, and then continuing 
N to the imaginary value —i/iR. This procedure is less well-defined, as it requires us to make a 
discrete parameter continuous. 
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3. Correlators from Matrix Quantum Mechanics 



3.1 Two-point functions 

We start by presenting formulae for the two-point function (T_ n / R T n / R ) to all orders 
in ^2 , from the Matrix Quantum Mechanics (MQM) approach. We will derive these 
formulae, valid at arbitary radius, starting from the infinite-radius formulae presented 
3l| . We start by quoting the closed- form expression for the infinite-radius two- 



ni 



point function (T_ g T g ), or more precisely the first derivative of the two-point function 
with respect to the cosmological constant, which is actually more convenient for our 
purposes: 



V r [i - W) r (f - ^ - Q) J 



(3-1) 



where q > 0. For clarity of presentation we will drop the leg-pole factors (T(—q)Y 
in what follows, keeping in mind that they can be restored whenever needed. 

Now we obtain the corresponding amplitudes at a finite radius R, using Eqs. ( |2.5| ' 
and (E13): 



(T- q T Q )n = R J— Im ( F 6 " * + <0 _ %t^±- 



sin 



where the first factor of R comes from the replacement of the ^-function by a Kro- 
necker 5 as in Eq. Q2.6I ). The differential operator in front is real and acts only on 
functions of /i, so it can be moved inside and we thus need to evaluate 

1 / T (I — ijj, + q) r (| — i//) 



2 sin (jgfyj y r (| - ifi) r (1 - ifi - q) 
This can be done very easily by expanding the operator as follows 



1 00 



- ) tj9,j 



j=0 

Using this we get the required expression as 




(3.2) 
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Next, we choose q = n/R. We see that the j th term from the first sum cancels the 
(j + n) th term from the second sum. So only the j — 0, 1, . . . , n — 1 terms from the 
second sum remain. Defining r = n — j, the above expression becomes 7 : 



(T_ n/R T n/R ) = Ree^ R J2^ ^ + ^ ^ 



-r(|-i M +(r-n-|)i) 



(3.3) 



In order to obtain the expansion of this expression in powers of l//i 2 , we can 
rewrite it in terms of the special functions: 



r(~ — i/i + a) r(| — iji — b) 



± 



(3.4) 



r(|-2/i + 6) r(|-z//-a) 

defined in Eq.(B.2) of Ref.|3J. We have: 

n/2 

(T_ n/R T n/R ) = Re e l ™/ 2fl ^ J* ((r - |)± (r - n - ±)±; //) , n even (3.5) 



r=l 



(n-l)/2 



Re 



e *™/«f|jH-(^,-^; M )+ E^ + ((r-|)i (r-n-|)i;/,)J, n odd 



Next we use the asymptotics for large /x: 

JF+(a, b- y) = e -^ a - b)/2 fi a - b f(a, b; pi) 



(3.6) 



where /(a, 6; //) is a power series in -j with real coefficients and starting with a 
constant term: 

/(a, 6; //) = 2 - i(a - 6) (a - b - 1) (3(a + bf - (a - b) - l) + O (J^J (3.7) 
It follows that, for even n: 

n/2 

(T- n/R T n/R ) = Re / ((r - (r - n - ^ 

r=l 

n/2 

= //^£/((r-|)i(r-n-±)l ;/ ,) 



r=l 



^_T(§-V*+(r- 



-r(|-z/i+(r-n-i)i; 



(3. 



7 Here and in what follows, we drop the R subscript in the correlators wherever it is obvious that 
they are at finite R. 
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The first step above follows because the function / is real. The final equality is 
true for all n, and not just even values. This then is the complete answer for the 
perturbative expansion of two-point functions of momentum correlators at arbitrary 
radius. 

Specialising to n — 1, we find the following expression, which will be useful later 

on: 

r (i - iii + 



(T-i/rTijr 



) 



After a T-duality 



R -> 1/R, fi-^/iR 



(3.9) 



(3.10) 



we get the unit-winding two-point function 



(T- R T R ) 



r(i-<Aii2 + f) 



r(i-^-f) 



(3.11) 



This expression was recently derived by Maldacena[19|]. We should note that the 
above answer has to be multiplied by the leg pole factor (r(— R)) 2 , which we dropped 
after Eq. Q. 



3.2 Four-point functions 



In this section we turn to the computation of higher point functions. In particular, 
we extend the results for the four-point function from MQM to finite R and then 
specialise to the case of unit winding modes. In this case we will be able to find an 
explicit all-orders result after summing an infinite series. 

Upto leg pole factors (which can be unambiguously restored when needed) the 
connected four-point function at infinite radius ispTf: 



d,{{T- q T q f) 



2 \ conn 

oo 



Im e 



'/Tvq 



r (2g — n + \ — iii) V (q — n + \ — ifi) 



r (i - z» 



r (-q + \ - ip) 



(3.12) 



where q > and the function T + is defined in Eq. (|3.4|). 
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Substituting Eq. (|3.4| ) in Eq. (|3.12|) we have 

r(|-z/i + 2g) r(|-z//) 



9,{{T- q T q Y) 



2 \ conn 



Im e 



ITTCj 



+ 



r (|-iAi + ff) 



+ 2 £ 



r (| - i/i) r (| - i/x - 2q) r (5 - i> - g) 
1)" /T(-g + 71) \ 2 /r (| - i/i + 2g - n) T(\-in + q-n) 



n=l 



n! 



r(-g) 



r (| - ifi) 



r (| - i/j, - q) 



(3.13) 



The connected finite- R amplitude is, therefore 



\2\conn _ D 2R M 



((T-qT q ) }ft nn = R ■ 



sm 



{{T-qT q ) 2 ) 



2 \ conn 
00 



We use the expansion Eq. (|3.2|) of the differential operator and set q = 1/Rto get 
{{T- l/R T 1/R y) co ™ = Re e*/* + ^ + 1 ^ ~ ^ + ^ 



F (|-^-2l?) r (l"^-2lj) 



^ i-D" m-i+n)Vr(|-2>+^-n) 



^ - ! V r(-i) ; r(±-<„-£) 



1(3.14) 



It is convenient to add and subtract a term corresponding to n = in the summation. 
This extends the sum from to 00, while the subtracted term changes the sign of 
the first term above, after which the first two terms combine into an J r+ . Thus we 
get: 



((T_ 1/R T 1/R fr nn = Re J*'* ( -i; A*) 



n=0 



-if /r(-i + n)\ 2 r(|-z At + 



nl 



r(-i; 



^ (2 ^ 2h) 



2i? J 1(3.15) 



The sum is now easy to evaluate using the integral representations for the three T- 
functions in the numerator that depend on n (see Appendix [B]). This finally leads 
to: 

2 N 



((T- 1/R T l/R ) 2 r™ = Re e^ R -fe fi) - 2 (±F 



1 77+/ 1 



(J- 

V2R' 



J_. 

2i?' 2W 
2 



/<) 



(3.16) 



One can verify that the two terms above are, respectively, the full (connected plus 
disconnected) correlator, and its disconnected part. 
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4. Correlators in the finite- N Normal Matrix Model 

Having obtained explicit expressions for all two-point and a particular four-point 
function from the MQM, as a function of the cosmological constant \x and radius R, 
we now attempt to recover the same results from the NMM. This first of all provides 
a test of the NMM and its effectiveness. But once we explore the systematics it will 
become clear that we can compute much more. In fact, we will obtain a complete 
combinatorial formula for the 2n-point functions of unit-momentum correlators. Via 
T-duality, this determines the corresponding winding correlators. We expect this 
to be useful in determining the full vortex condensate to all orders in perturbation 
theory. 

As mentioned before, in the process of studying the NMM we will encounter a 
rather surprising result: for the purpose of computing correlators, one can actually 
take N to be a small finite value and yet obtain the correct answer as a function of 
fi. The finite value of N will be determined by the operators whose correlators we 
are calculating. For this purpose it is convenient to classify tachyon correlators into 
sectors labelled by an integer, the total positive momentum P flowing through that 
correlator, measured in units of 1/R. For example in (T^k 1 /RT-k 2 /RT mi /RT m2 /R) , 
where ki, k 2 , mi, m 2 are all positive, the total positive momentum is P = mi + m 2 = 
ki + k 2 . This number will determine the minimum value of N required in the NMM 
to compute these correlators. In what follows we will first consider all correlators in 
the sectors P = 1 and P = 2. In the former case there is only a single two-point 
function, while in the latter case we have two, three and four-point functions. After 
presenting some examples we will discuss why the theory works in this way. 

4.1 Two-point functions: examples 

Example: n = 1 

We begin by computing the two point function of the unit momentum operator. 
Since total momentum is conserved, this operator is paired with the one of negative 
unit momentum. So we will calculate the two point function (T_i/ R Ti/ R ) of unit 
momentum operators. 

We first calculate the partition function of NMM at N = 1: 




(4.1) 
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Setting 2 = y/me l6 i dzdz — > dmdO, we have: 



poo p2tt 

Z%M l M {t = 0) = / / dm d6 e — io, 
Jo Jo 



og m 



/■oo 

2tt / rfmm(^ 1+ ^)e-^ 



o 

= |,-(^--)r(, + |-^) (4.2) 

As a function of u, this is not the correct partition function of the c = 1 string, but 
it reduces to the correct partition function if in the above expression we set v — 
and compare this with Z c= i(j^, t = 0,t = 0). This fact is a direct consequence of the 
claim in Ref . pi | , see Eq. ( p. 15 ). It is also worth noting that the partition function 



at iV = 1 is not invariant under T-duality. In fact, T-duality in the NMM partition 
function is recovered only in the limit iV — > oo. This makes it clear that the correct 
partition function, as a function of fi and R, can never be recovered at finite N. 

For correlators, things are quite different, as we will now see. For the two-point 
function, we find: 

POO p2tT 

rN=l (+ _n^_ / J^jn^ c -vrnR+(Rv-l+^)\o gm 



poo pzir 

d-idiZ^^l M {t = 0) = / / dmdOme' 

Jo Jo 

= |,-(-^)r(, + i + ^) (4.3) 



From Eq. ((O) and Eq. (O) we have: 



in^/ M (t = o) = »H ^fr + ; + y) (4.4) 

r ( Z/ +2- 2R) 

Finally, we have to analytically continue v = —ifi. The result is complex, but can 
easily be seen to have the form of an overall phase times a real power series in l//i 2 . 
Dropping the phase is then equivalent to taking the modulus of the above expression. 
This gives: 



(T-i/rTi/r)nmm - V R 



1 (I - ^ + 27?) 



(4.5) 



r — ifi — jr) 

which agrees with Eq. (|3.9|) upto the prefactor, which indicates that the 

"tachyons" of the NMM are normalised differently from those of MQM. Indeed we 
will argue later that the relationship is: 

T n /R\NMM = V> n ^ R T n /Ft\MQM (4.6) 
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We have discovered the surprising result that the exact two-point correlator of 
unit momentum tachyons is correctly calculated (as a function of fi and R) using only 
the lxl Normal Matrix Model! According to Eq. ( |2.15| ), we should have expected 
the result to be correct only for fi = i/R. We will see that a similar feature holds for 
all two-point correlators, though the minimum required value of N depends on the 
correlator under consideration. Later we will extend this observation to higher-point 
correlators. 

Example: n = 2 

We consider another example, the correlator (T_ 2 /_r72/k)- in this case, according 
to the prediction in Eq. ( |2.15|) , we can perform a calculation at iV = 1 and the result 
so obtained will be valid at the special value of the cosmological constant v — 1. 
However, we now face a puzzle. In the NMM at N = 1, one cannot distinguish the 
four correlators: 



(T-2/rT 2 /r), (T-2/rTi/rTi/r), (T-i/rT_ 1 / r T 2 /r), (T-i/rT^i/rTi/rTx/r) 

(4.7) 

because all of these are represented by the same NMM correlator (z 2 z 2 ). Therefore, 
assuming Eq. ( |2.15| ) continues to hold, either it has to be the case that all four 



correlators become the same at v — j^, or else at best we can only hope to obtain 
some linear combination of them. 



The calculation is straightforward and upon continuing to v 
the modulus, we find: 



rp \N=l 

2/RJ-2/R)nMM — V 



-2/R 



r(I 



r (2 ^ 2r) 



This can be compared with the known result from Eq. 
present case, and changing to the NMM normalisation via Eq. 



-i/j, and taking 
(4.8) 

Specialising to the 
gives us: 



(T-2/R.T2/R 



) 



r( ifi + „ 



+ 



r( 2 i^i 2r) r( 2 %fi 2R/ 



(4.9) 



Comparing Eqs. (|4.8|) ,( [4.9|) , we see that the NMM result for this correlator at iV = 1 
is not correct. This is not a surprise. But now we see that it is incorrect even at 
the special value fi = i/R, which appears to contradict Eq. fl2.15|) . As we will see, 
this is due to the fact that the same NMM correlator can describe different tachyon 
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correlation functions for low N. Indeed, one can check that the answer we have 
obtained at N = 1 in Eq. ( |4.8| ) is actually a linear combination of the correlators in 
Eq. ( |4.7|) as calculated from matrix quantum mechanics. 

Let us continue by evaluating the NMM correlator at iV = 2. In this case 
the operator we are dealing with is T2/R ~ tiZ 2 which is linearly independent of 
(Ti/r) 2 ~ (trZ) 2 once Z is a 2 x 2 matrix, so there is no longer a risk of mixing for 
the operators in Eq. (|4.7|) . The computation is given in an Appendix, and leads to 
the answer Eq. (|A.5|) , which after changing to the NMM normalisation is: 



(T~2/rT2/r)nMM 



r(|-<A* + s) + r (i-v + s) 



r (2 2i?) I" (2 i[l 2 r) 



(4.10) 



Following Eq. ( 2.1 5| ) we would expect that this should give the correct answer for 



H = 2i/R. But now there is a surprise, since in fact it agrees perfectly with the MQM 
result Eq. ( |4.9|) for all values of li. Thus for the purposes of calculating (T_ 2 /ijT 2 //j) 
in c = 1 string theory, to all orders in the string coupling, a 2 x 2 matrix model is 
sufficient. 

To summarise, we have found evidence that an NMM calculation of tachyon 
correlators at finite N (where the minimum required value of iV depends on the 
correlator in question) gives the correct tachyon correlators for the c = 1 string, to all 
orders in perturbation theory. Below we will collect more evidence for this property, 
which appears to go far beyond the result of Ref. [21] as stated in Eq. ([2.15|) above. 



4.2 Two-point functions: general case 

Let us now consider the general case (T_ n / R T n /R) . and try to derive this result from 
the NMM. We will find that for this correlator, the NMM with N = n is sufficient 
to give the correct result. Indeed, when we compute in the N x N NMM starting at 
N — 1 and increasing N in integer steps, we obtain the right c = 1 string correlator 
(as a function of li) as long as N > n, though not for N < n. Thus the NMM 
calculation "stabilises" at a certain minimum value of N. 

Since we will be computing normalised correlators, we start by computing the 
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(unperturbed) partition function at a general value of N. This is given by 



/»oo 1 ■ rzn 1 v 

0) = / Y[dm r / JJd0 r (4.11) 



AT 



j<fc 

The next step is to perform the integration over the #'s. In general this will be quite 
tedious, because one has to pick out terms which are independent of 6 by expanding 
out the Vandermonde factor. However, we notice that since the above expression 
is invariant under permutations of m's, we can determine all terms surviving the 9 
integrals if we know just one of them, by permuting the m's among themselves. 

The first such term is just the product of the first term from each of the Van- 
dermonde factors, which is m^ -1 ??^ -2 • • • itin-i- Thus we have, after evaluating the 
9 integrals 

„oo N N-l 

Z N NMu{t = 0) = (2k) n N\ / J] dm r J] mf ^ (4.12) 

J ° r=l j=l 

x e -KS^ =1 <) + (^-iV + « f i)(E^ 1 logm r ) 
N 

= (2n) N N\ ft v-(»+Hr-mr(v + I - ( r - I) l) 

r=l 

^From now on we will restrict to the case N = n. 

The next step is to compute the two point function and then normalise by the 
above partition function. We have: 



z i 



/ib ib iiv \ I 

n^*ni*i-**i 2 £*r £ 

r=l j<k \l=l J \l=l 

Xe -"(Z? = l(ZrZr) R )+(Rv-n+K=±)(Z: =1 logZ r Z r ) 
poo n p2tt n n 

= ] J dm r / II d6 r Y[ {mj + m k - y /m j m k (e ie ' k + e~^' fc )) 

Jo r=l J ° r=l j<k 

n \ / n 

E(^) n ^ £(v^,) n e- 

r=l / \r=l 

Xe -"(E"=im?)+(^-«+¥)(S = i tag"*) 

In this case also we can avoid tedious calculation by applying the permutation 
trick. The contribution to the first term from the Vandermonde is same as be- 
fore, and the contribution from tiZ n tiZ^ n is ^™=i m r- The net contribution is then 
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(Sr=i m r) m i lm 2 2 ■ ■ ■ m n _i. Proceeding as before we have after the 9 integrals 

/oo n / " \ n—l 

l[dm r lj2m n r ] I] 



n-l 



r=l \r=l / j=l 
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r=l 



(4.14) 



From Eq. (|4.12|) and Eq. Q4.14 ) we find (after changing variables j — > n + 1 



(t.^t^w _ , ^ r( i_ i/1 + (r _ n _i ) i ) ( 4 - 15 ) 



As before, we analytically continue z/ = —ifi and take the modulus to get: 

r(|-^+(r-i)i) 



irp rp \N=n —n/R 

\ 1 -n/R 1 n/R/NMM ~ P 



E 



(4.16) 



-r(i-v+(r-n-l)i) 

After changing normalisation via Eq. ( f4.6| ), we see that this agrees perfectly with 
Eq. (pD. 

The above calculation was performed with matrices of rank N = n. It can 
easily be repeated for the other cases. When N is smaller than n, we find that the 
answer, as a function of /z, is not equal to the correct two-point function, and does 
not become the correct one even after choosing /j, = in/R. As before, this is due 
to "contamination" by correlators of higher point functions carrying the same total 
momentum, because for N < n the corresponding correlators in the NMM are not 
all linearly independent. For N > n, instead, we actually get the same final answer 
as for N = n. The calculational procedure we described above seems to suggest that 
extra terms arise for N > n, but actually they are cancelled by contributions from 
the 9 dependent terms in the Vandermonde factor. Thus when we take the ratio of 
d^ n d n Z and Z we end up with the RHS of Eq. ( [4.151 ). Therefore as long as we take 
N > n, we get the right answer (independent of N) for every N. This is what we 
referred to as "stabilisation" above. 

4.3 Four-point functions 

Now we would like to compute the four-point function in the Normal Matrix Model. 
For N—l, the calculation has already been performed, since as we noted above, it 
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is the same as the corresponding calculation for the two-point function in Eq. (|4.S| ) 
(more precisely the disconnected four-point function is the same as this two-point 
function). As we explained there, the result so obtained is a linear combination of 
the correct two, three and four-point functions of the c = 1 string, and to distinguish 
them we need to go to a higher value of N. Accordingly we have computed the above 
four-point function using the N = 2 NMM. The derivation can be found in Appendix 
[C| and the result is: 



({ T -1/rTx/r) 2 )nMM — fx 2 



T+(Z L. u ) - ±(t+(A L.»yy 

J \2W 2R ' r 1 ) 2y V 2R ' 2R 1 W J 



(4.17) 



Changing from MQM to NMM normalisation using Eq. (|4.q) , and inserting the usual 
1/R factor, we see that Eq. (|4.17|) above is identical to Eq. ( |3.16| ). 

For completeness, let us briefly consider the two three-point functions 

(T-2rTrTr)nMMi (%rT-rX-r) NMM ( 4 - 18 ) 

The two are actually equal to each other because of the symmetry X — > —X, where 
X is the Euclidean time direction. We have calculated these correlators both from 
MQM and NMM (at N = 2) and the agreement is exactly as for the cases considered 
above. 

4.4 Why it works 

As we reviewed in Section ^|, the Normal Matrix Model determines every momentum 
correlator by differentiation with respect to the momentum couplings t, t. However, 
the correlators so obtained should only be correct in the limit N — > oo ("Model 
I", Eq. fl2Tp ) or the special values N = vR ("Model II", Eq. (|2~T5|)). Now in the 
previous subsections we have shown in several examples (including the infinite set of 
two-point functions) that, given the total momentum P flowing in the correlator, the 
NMM with matrices of any rank N > P suffices to compute the correlator completely 
as a function of fi and R. In view of this, the NMM appears to go beyond its expected 
range of validity. Here we will give an explanation as to how this comes about. 

The basic observation is that the phenomenon we are observing is not to be 
viewed as an application of Model II, but rather of Model I. Indeed, using Model II 
and a definite value of N, it is clear from Eq. (|2.15|) that the answers obtained are 
correct only for a definite value of u, namely v = N/ R. This relation between N and 
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v defines a line in (N, u) space, and the points on this line where N takes integer 
values are the ones where the procedure works. However, it is clear that in this way 
one can never recover the full v dependence at a fixed N. 

In contrast, in Model I one is supposed to compute correlators at an arbitrarily 
large value of iV and in the limit N — > oo, the correct answers are obtained as a 
function of v. What we will now show is that, after computing a given correlator of 
total momentum P in this way, and then dividing by the partition function, infinitely 
many terms cancel out exactly in the ratio. The remaining terms, which actually 
contribute to the correlator of interest, are the same as one would compute for a 
finite value of N, namely N = P. 

The argument goes as follows. From the derivation we have given in the previous 
subsections and the appendices, any correlator is generated (after 9i integrations) 
by inserting an expression of the form Yl%=i m T va ^° ^ ne m « integrals, where 
correspond to ordered partitions of P. Therefore we should first of all choose iV large 
enough so that all such partitions can be realised and are distinguishable. This is 
possible for N > P. For iV < P we will miss some partitions, and thus the answer 
cannot be correct. But the case N > P realises the same partitions as the case N = P 
and thus gives the same answer. This causes what we earlier called "stabilisation", 
which amounts to saying that the result for iV = P is identical to the result for any 
N > P, and therefore for N = oo. Invoking the converse of stabilisation, we can 
therefore start with the model defined at N = oo and "bring back" the value of N 
to any finite value N > P without changing the result. This explains why a finite- N 
matrix model is sufficient to compute any momentum correlator. 

4.5 Combinatorial result for 2n-point functions 

We have shown that the NMM is an effective tool by re-computing known correlators. 
Now that we understand how and why it works, we apply it to compute a new result: 
the full (connected plus disconnected) 2n-point function ((T_ii R Ty R ) n ) for every n 
and to all orders in perturbation theory. The result, derived in Appendix [D|, is the 
following: 



{ki} i=l 



n 



r(i-v-(i-s)s) 



(4.19) 



21 



with C({ki}) defined as: 

^})=B-if n 



i=i 



fa - Pi 



(4.20) 



Here, {k{\ are strictly ordered partitions of n(n + l)/2, namely: 

n(n + 1) 



h > k 2 > ■ ■ ■ > k n , 



(4.21; 



i=l 



and denote permutations of the n numbers n — 1, n — 2, • • • ,0. 

Let us examine this result more closely. In principle, for every n the answer is 
a sum of terms, each one being the ratio of n T-functions divided by n T-functions. 
However in practice, some of the numerator and denominator terms can cancel out. 
We can see this more explicitly if we list the first few special cases, of which the first 
two have already been noted above: 



(T-i/rTx/r) 
((T- 1/R T 1/R f) 
{{T_ l/R T l/R f) 



T (2 2r) 
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(4.22) 



W + 2R' 



+ 4 



T{\-i^ + ^ R ) r(|-z//+^ ; 

r(i 



r( 2 %\i 2R ) r( 2 ifi 2R ) r( 2 %n 2R) 

The pattern emerging so far is misleadingly simple, as we see with the next example, 
the 8-point function: 
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2/i 2jR/ 
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(4.23) 



We see that as the number of operators in the correlator grows, one gets products of 
more and more T-functions in the numerator and denominator. In this example we 
also see clearly that the coefficients are perfect squares. 

Ideally one would like to know the connected part of the 2n-point function. 
In principle this can of course be obtained by repeated application of Eq. ( |4.19| ), 
but one would like a more explicit and useful expression. However, for the most 
likely application, to the vortex condensate, we will not really need to make the 
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distinction between connected and disconnected correlators. The vortex condensate 
corresponds to the partition function of a perturbed theory, and to find the connected 
component of that it suffices to take a logarithm. We will discuss this issue further 
in the following section. 



5. Applications 

5.1 T-duality at c = 1 

In this subsection we discuss how our results can be applied to check T-duality of 
the c = 1 matrix model. As we have seen, in the Euclidean (finite-temperature) 
MQM, the momentum and winding modes with respect to the time direction are 
independently defined. The former arise from macroscopic loops defined in terms of 
fermion bilinears, while the latter are Wilson-Polyakov loops in the thermal direction, 
which project the theory onto nonsinglet sectors. From the continuum description 
we expect that there should be T-duality between these two sets of observables. 
Indeed, in Ref.||16|| it has been formally argued that, like the momentum-perturbed 
matrix model, the winding-perturbed MQM also corresponds to the r-function of a 
Toda hierarchy. To understand T-duality better, one would like to compare explicit 
correlation functions computed from the momentum and winding sides. 

An attempt to directly check T-duality was made by Maldacena in |19|], where 
the following two quantities were compared: (i) the two-point function of unit- 
momentum tachyons, after T-duality, and (ii) the partition function of MQM in 
the adjoint sector. From Eq. (13 . 1 If) we see that (i) is equal to: 



(T- R T R ) 



r(i-i/*fl + f) 



r(|-^-f) 



(5.1) 



However, at this point we recall that leg-pole factors of r(— \q\) were dropped after 
Eq. (|3.1| ). Restoring them and taking the large-// asymptotics of this correlator, we 
find 8 : 



(T- R T R ) = (T(-R)) 2 ^R) R U + ^(R- ^)^ 2 + O O^- 4 )) 



(5.2) 



On the other hand, (ii) is obtained by solving MQM in the adjoint sector. In the large 
iV limit, Maldacena obtained the leading (tree level) contribution to the partition 



i The factor R R was not written in Rcf.|19|. 
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function in this sector as: 



- (Wadj) = — ^—n R = (r(i?+ i)r(-J2)) V (5-3) 



Z sing x 48^71-^ 4vr 

The power of /i agrees with that in the leading term of Eq. (|5.2|) . The remaining dis- 
crepancy can be assigned to the normalisation of the fundamental Wilson-Polyakov 
loop (or equivalently to the normalisation of the original momentum modes), and we 
see that Eqs. Q5.3p and ( |5.2j ) agree to leading order if we change the normalisation of 
this loop variable to: 

This is a relatively simple change of normalisation 9 , and appears to specify the basis 
in which T-duality holds in MQM. 

It is not entirely surprising that one needs to change normalisation of the matrix 
model observables in order to implement T-duality. Indeed, this duality is most man- 
ifest in the worldsheet or Liouville approach, in which the momentum and winding 
vertex operators come with a natural normalisation and are related to each other 
by the simple change (X L ,X R ) — > (X L , —Xr). On the matrix model side, momen- 
tum operators in the MQM are related to the corresponding Liouville operators by a 
change of normalisation, Eq. ( |2.4| ). So one should expect that winding operators in 
MQM are also related to Liouville winding operators by a change of normalisation. 

This is not to say we understand the nature of these normalisation factors in 
general. In fact, as stressed in Ref . 1 19j , we need more examples in order to check 



the consistency of this picture. As an example, if one could compute the genus-1 
correction to the adjoint sector partition function, this could be compared with the 
genus-1 term in Eq. ( |5.2j ). Similarly, if one could compute the leading term for those 
higher representations that correspond to 2n-point functions of the winding tachyon, 
then one could match this with the asymptotics of the latter, which can be read off 



from our results in Section 4.5 



There will also be representations corresponding to the correlators of multiply 
wound tachyons T nR . These correlators can be found by T-dualising the relevant 
momentum correlators, for example the two-point functions are found by T-dualising 



9 Notice that the normalisation factor becomes trivial at the special radius R = 1. 
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Eq. (TO , leading to: 
(71 ni? T n/? ) = (T(-nR)) 2 

, \ ft / ^2/ nR(nR -l)((n 2 - l)R 2 -nR-l) „ , ,,\ 
= n(/^r*(r(-ni*)) 2 ( 1 ^ Al 24fi2 J A* + ^ (^ 4 ) J 

In the matrix model, this should correspond to the fundamental Wilson-Polyakov 
loop with a contour that winds n times over the time direction. In principle we are 
allowed an independent choice of normalisation for each winding number. In fact 
the momentum and winding modes have corresponding freedoms in normalisation, 
and the only thing relevant for T-duality is the relative normalisation between them. 
So when we consider the nonsinglet sector related to multiply wound loops, and the 
corresponding tachyons of n units of momentum, the leading-order comparison will 
be used to fix the normalisation and the loop corrections will constitute a genuine 
check of T-duality. 

To summarise, we have not been able to address the problem of T-duality but 
only set up one side of it. Namely, we have exhibited the all-orders finite-radius 
correlators computed from the momentum side, after performing a T-duality trans- 
formation. This constitutes a prediction to be checked once it is properly understood 
how to perform nonsinglet computations for different representations and to higher 
orders in string perturbation theory. 

There is one more intriguing point that we would like to mention. The correlators 
we have computed take very special values at the self dual radius R = 1, the point 
of enhanced SU(2) symmetry. In particular, all loop corrections to the two-point 
function of unit momentum tachyons vanish, as can be seen from Eq. ( |3.11|) . Thus 
the tree level answer is exact 10 . By T-duality the same property should hold for 
the two-point function of unit winding modes. It is plausible that one could extract 
this simple property just from the structure of the nonsinglet Hamiltonian - in this 
case it is the adjoint Hamiltonian that was studied in Ref.|19|, specialised to R = 1. 
Similarly, at R = 1 the other two-point functions have perturbation series that 
terminate at a finite number of loops, as one can easily check from Eq. ( |5.5|) . So, 
for consistency this must also be a property of the antisymmetric-antisymmetric 
10 This was already known long ago, for example as the puncture equation in the Kontsevich- 
Penner model ]22|. 
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representations referred to above. It may be simpler to derive this kind of general 
result in the nonsinglet sector than to actually compute coefficients with precision. 



5.2 Vortex condensate and black holes 

It is believed that the Euclidean 2D black hole background, defined in the continuum 
by an SL(2, R/U(l) CFT, is equivalent to the c = 1 matrix model perturbed by 
fundamental Wilson-Polyakov loops: 

Smqm -> S M qm + A W N + A >% (5.6) 
The basis for this belief is the FZZ conjecture! Hl> which relates the black hole back- 



ground to Sine-Liouville theory 11 . Via the equivalence in Eq. (|2.8|) , the latter is the 
same as the perturbed background above. 

To be precise, the FZZ conjecture is not really an either/or statement wherein 
one uses either the black hole background or the Sine-Liouville perturbation. It has 
increasingly become clear that the backgrounds that one might call "black hole" or 
"Sine-Liouville" are the same, and both perturbations are turned on simultaneously. 
Depending on the value of the worldsheet coupling, one or the other of these per- 
turbations is more dominant, but for example the exact correlation functions have 
poles corresponding to both perturbations 12 . In the present work we will not focus 
on these details, but will be content to treat the black hole story as a motivation to 
understand the vortex condensate: 

/ e XW N +XWff\ 

MQM 

One way to compute this condensate would be to sum over an infinite set of nonsinglet 
sectors in the MQM with some definite weights. However, as we have seen, the 
technology to do this seems rather limited at present. An alternative is to assume 
T-duality to compute the correlator: 

OO OO \ n \ 771 I \ I QiTl 

<e AT * +xr -*> = E E ^—^rT^-rD = E rU(( r « r ^) n ) ( 5 - 8 ) 

z — ' z — ' n\ ml z — ' (n!) 

n=0 m=0 n=0 v ; 



11 This conjecture has been proved by Hori and Kapustin[|33[ in the J\f — 2 supersymmetric case. 
As Maldacena has argued [^9], suitably orbifolding both sides of their argument leads to a proof for 
the bosonic case. 

12 See for example Ref. Q. We are grateful to Ari Pakman for explaining this to us. 
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where the last equality follows from conservation of winding number. 

Now from the computation in Appendix [B], we have the following result after 
T-duality: 



{ki} 



-T r (~ - ijiR 



(5.9) 



where {ki} are strictly ordered partitions of n(n+ l)/2, and C({ki}) are the combi- 
natorial coefficients given in Eq. ( |L).9| ). 

The above correlators contain both connected and disconnected contributions. 
We can now pass to the generating function: 



(e 



XTr+XT_ 



'> = £ 



I A 



2n 



n=0 



in 



\)2 



i=i 



k - \)R) 



r (\-i^R-(i-\)R) 



(5.10) 



This is the partition function in the presence of a vortex condensate, and its logarithm 
is the free energy of the perturbed theory. So one does not need at any point to 
compute individual connected correlators. 

The above expression is completely explicit and does not require integrating any 
equation or developing a recursion relation. We expect it will be useful to to extract 
physical quantities of interest related to the Euclidean 2d black hole. This is beyond 
the scope of the present work, however, and we hope to return to a more detailed 
analysis of this formula in the future. 

Again it is worth pointing out that at the selfdual radius R = 1 the vortex 
condensate is known exactly, though deriving it from the above expression would 
not be the easiest way. The puncture equation of Ref. f22|| simply tells us that: 



(e AT *+ AT ->| R=1 



-ifiXX I 



(5.11) 



and one can check easily that this agrees with the cases in Eq. ( |4.22| ) specialised to 
R = l. 

The significance for the Euclidean 2d black hole of this simple result has not, to 
our knowledge, been explored. While it is true that the black hole CFT corresponds 
to a radius R — |, it is believed [03] to have a marginal deformation to other radii 
at least in the range 1 < R < 2. So the physical consequences of the simple formula 
above at R — 1 would be worth understanding better. 
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6. Conclusions 



In this work we have examined the familiar c = 1 bosonic noncritical string theory, or 
rather its Euclidean (finite temperature) version, from the perspective of correlation 
functions. Both old and new techniques were used to develop simple, elegant and 
explicit formulae as functions of two variables: the cosmological constant fi and the 
compactification radius or inverse temperature R. The key results are summarised 
in Eqs.( [3.8|) ,( 3.16|) , (|4.19| ). In addition we have shown that the Normal Matrix Model 
is a powerful computational tool. 

An obvious extension of this work would be to the case of noncritical type 



strings [ 34], [35 1 . In Ref. [36], explicit expressions are obtained for the partition func- 
tions of type OA and OB strings in the presence of fluxes. These expressions are richer 
than the corresponding ones for the bosonic noncritical string, both because of the 
flux dependence and because they are nonperturbative in /i. Our work should gen- 
eralise quite straightforwardly, particularly to the Euclidean type OB case, and the 
correlators so obtained will contain nonperturbative information about the theory. 

A detailed investigation into the physical questions that motivated the present 
exercise, namely a better understanding of the 2d black hole background as well as 
of T-duality in the matrix model, is left for subsequent work. We also note that the 
physical origin of the Normal Matrix Model has not yet been understood. As it is 
clearly a correct and useful description of the c = 1 string, and moreover makes sense 
only in the Euclidean context, it would be worth trying to put it on a similar footing 
as MQM in terms of the dynamics of some appropriate (Euclidean) D-branes. 
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A. Computation of two-point functions in the NMM 

Here we present some of the details of how to compute two-point functions in the 
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Normal Matrix Model. To start with, for the partition function we have 

^NhlM^ = 0) = J d 2 Z\d 2 Z2 \Z\ — Z 2 \ 2 

As before, we change variables Z{ = ^JrrTit 1 ® 1 , d 2 z; L — > drriidOi and we get 

Znmm^ = 0) = / dm x dm 2 / dd x dd 2 (mi + m 2 - v/mTmi^ 12 + e^ 12 )) 
Jo Jo 

xe -u(m R +m R )+(Ru-2+ ^p- ) ( log mi + log m 2 ) 

poo 

= Air 2 / dm x dm 2 [m x + m 2 ) (m im2 ) (^-2+^f i ) e -K+0 

= S--^ + --)r (* + 1 - 3L) x ^-H-^r („ + 1 - x) , (a.2) 

where #12 = 61 — 62- In a similar manner we have 

d^ 2 d 2 Z^ 2 M (t = 0) = y rf 2 ^rf 2 z 2 |z a - z 2 | 2 (^ 2 + 2 2 )(% 2 + ) 

xg -i/((2;i5i) fl +(22S2) B ) + (fii'-2+^-)(logzizi+log22Z2) 

/*OC /*27T 

= / dm 1 dm 2 / c^cf^ (mi + m 2 - v /m 1 m 2 (e^ 12 + e^ 12 )) 
Jo Jo 

x ( m 2 + m 2 + mim2e 2^ 12 + mim2e -2^ 12 ) e -, (m f +m «) + (i ? ,-2 + S f i)(log mi + lo gm2 ) 

/•OO 

= 4vr 2 / c/m lC /m 2 ( mi + m 2 )(m 2 + m 2 2 )(m x m 2 )( Ru - 2+E ^) e~ v ( m * +m *) 
Evaluating the integrals on m x , m 2 we get 

d_ 2 d 2 z^ 2 M (t = o) = g„-M + &)r (v+l + A) x i/-W-*)r (* + 1 - 

+ S^ (H_i+A) r (^ + l + 2i?) x ^-*)r (1/ + 1 - i) 

(A.3) 



From Eq. ( A.2 ) and Eq. ( A.3 ) we have 

(A.4) 



/T rp )N =2 : .-2,r r{v + \ + TR) 1 r ^+I + ^) 



As before, to get the correct two point function we have to analytically continue 
v = —ifj, and take the modulus of the above expression. This gives: 



/rp rp \N=2 ,,-2/R 
{-L-2/R-L2/R/NMM - V 



r(|-*> + ^) , r(|-i/* + ^) 



^ (2 ^ 2r) r ( 2 «/i 2 ^) 



(A.5) 
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B. Evaluation of a summation in the MQM four-point func- 
tion 



In order to show the equivalence between Eqs. ( 3.15 ) and ( |3.16| ) we need to prove 
the following identity: 



£ 

n=0 



r?! 



r (-i) / r(i-i M -i) lr(i-^-^) 



(B.i) 



Let us start with the expression: 



1)" /T(-B + 



r(-i 



r (| - 1> + a _ „) 



(B.2) 



n=0 "" v * v R' 

Using the integral representation of the V function we write this as: 



£ 



The sum over n can now be performed immediately and we have 



(B-3) 
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dH e *3 (i^) - * ~ x tJ 



(B.4) 



Using the change of variables t x — > + I 2 -) and performing the integral on t\ we 
get: 

£ = — ^- [ dH tj* -1 tl'^^' 1 (t 2 + t 3 )* e"* 2 "' 3 (B.5) 



n-h > 

We next introduce a parameter a which allows us to write the above equation as 



da 



r<-i 

Changing variables ti — *> at{ we have: 
1 



(B.6) 
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a=l 



Using the relation: 



d_ 

da 



a 
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r(— n + m) 

n-n) 



(B.7) 



(B.8) 
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we finally have: 



£ 



(r(|-vi + ^)V 



%H 2R ) 

Using Eq. ( |B.2| ) and dividing both sides by T(| — i\i 
Eq. Q. 



(B.9) 

we immediately get 



C. Four-point function in NMM 

We now briefly describe the calculation of the connected four-point function of unit 
momentum modes in the NMM. This is obtained by differentiating the free energy 
T with respect to the couplings. We thus have: 



((T- 1/R T 1/R f) = d\d\T 

= ((T_ 1/J? T 1/fl ) 2 ) di ™ - 2(T_ 1/R T 1/R ) 2 



(c.r 



where T = In Z^mm- The second term in the above equation can be calculated from 
the NMM with N = 2 and is given by: 



((T- 1/R T 1/R y) 



2 \ disconn 



<(trZt) 2 (trZ) 2 > 



N=2 
NMM 



V 



-2/R 



The explicit calculation is very similar to the calculation of (T-2/rT2/r) from the 
NMM. The disconnected piece is simply the square of the two-point function listed 
in Eq. ( |4.4j ). Putting everything together the connected four-point function is given 
by: 



((T-i/rTi/r) 2 )nmm - v 



-2/R 



-2 



(C.2) 



Analytically continuing v = —ifi and taking the modulus, and then using the 
definition of T + in Eq. (|3.4j), we finally get: 



{{T-i/rTi/r) 2 Tnmm = (aO 2/R 



T+( JL L. u ) _ J L- w ) 

^ V2i?' 2Ji'W 2^ V2fi' 2i?'W 



(C.3) 
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D. 2n-point functions in NMM 

Here we present the detailed calculation of the 2n-point functions from the NMM. 
In what follows we will take the rank of the matrix, N, to be equal to n. We have: 



I Zi Zj | 



(d^rz^t = o) = I n^n 

J i=l i<j 



|2 

i=l i<j \ i=l / \i=l 




We would now like to make the substitution z% = y/rfii e l6i and perform the 6 integrals. 
The remaining integrand will then be a function of the rrii and we will find that it has 
the form (^2{ k .y C({hj}) 2 f| f + permutations j e _5jVMM . Here {ki} are positive 
integers corresponding to strictly ordered partitions of n(n + l)/2, i.e.: 

n 

J2k l = r -^, k 1 >k 2 >--->k n >0 (D.2) 
i=i 

The permutations referred to are of the mj. Because the are integration variables, 
summing over permutations simply amounts to multiplying by a factor of n\. The 
constant coefficients have been labelled C({/c;}) 2 in anticipation of the fact that they 
will turn out to be squares. After performing the integration over m 8 and dividing by 
Znmm we get the final answer as a sum of ratios of products of gamma functions, with 
each term in the sum corresponding to a strictly ordered partition {ki} of n(n + l)/2. 

We will first show that the coefficients are perfect squares C({ki}) 2 . After that 
we will turn to the calculation of the C({ki}). Consider the expression: 




The full integrand is then UU times the exponential factor. Because the action is 
independent of the 0's, the entire ^-dependence of the integrand is in UU. Note 
that U has only positive powers of e ldi and U has only negative powers. Only the 
^-independent terms in the expansion of UU will survive the 6 integrals. 
It is easy to see that if we expand U,U then we get: 

n 

U = Zi 1 + permutations 

K} i=i 



n 



U = Y^ c (i k i}) JJ z k i l + permutations, (D.4) 

K} »=i 
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with {ki} defined as before. It is now clear that the coefficients of ^-independent 
terms in TT must be perfect squares, as the phase of a term in the first expression 
of Eq. (p. 4\ ) can only be cancelled by the complex conjugate term from the second 
expression, which has the same coefficient as the first term. 

Let us now determine the coefficients C({cti}). First we note the following prop- 
erty of the positive phase part of the Vandermonde: 



a, 



j<k V j=l 

where V is a particular permutation of the n integers {n — 1, n — 2, ■ ■ ■ ,0) and Vj 
denotes the j th element of the permutation V 13 . The sign for the first permutation 
is positive by construction. Any other permutation can be arrived at by a series 
of interchanges Zi <-» Zj. Each such interchange introduces a minus sign in the 
Vandermonde. Thus even permutations have a positive sign, while odd permutations 
have a negative sign, leading to Eq. ( p.5|) . Expanding the first factor in Eq. (|D.3| ) 



in a multinomial series and using Eq. (p.5|) we get: 




meiii r^K'iiB-i) p n 

V j=l 



n 

V 



{ft} V 4=1 V Pi J 

where are the unordered partitions of n. 

Let us examine the possible values of the exponent ki = $ + V{ in the above. If 
ki = kj for some i ^ j then the corresponding coefficient is zero. This can be traced 
back to the fact that the expression Eq. (p.3| ) is odd under pairwise interchange of 



the z's. Therefore we can rewrite the above as: 

u= e E(- l ) p ^^' ! " E ^ < ^~ p^) )^ ^ '• (D7) 

ki^kj V 4=1 \ Ki Hi J i=1 

E t fc,=n(n+l)/2 

Because the fc, are all distinct, we can limit ourselves to strictly ordered sets satisfying 
k\ > &2 > ■ ■ • k n . The other orderings are obtained by permuting these ones, or 



13 



For example, Vj = n — j when V is the identity permutation. 
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equivalently by permuting the Zj's. Thus we have: 

n 

U = ^ C ({ki}) TT z- 1 + (permutations of zi) (D.8) 

fci>fc 2 >--->fc„ «=1 
£.fc i= n(n+l)/2 

with 

P i=l \ fa Pi J 

Finally we combine U with U and integrate over the angles to get: 



/n n 
a — 1 n. l a — i 



)»■*>■:".. = lW I TT ,/„, Vrili ..lV-TT m * ! ' e 5:r=i(-'""f+(i&/-n+%i)logm i ) 

i=l i=l 

+permutations (D.10) 

n 

(2n) n n\ C{{h}f J] zy-(5+^-™+l)^)r (I + ^ + - n + 

4=1 



Using the expression for the partition function Z^mm from Eq. (|4.12|) for N = n we 
have: 

(d_ 1 d 1 ) n Z^ I n M _ n / R ST r(Sh " 2 TT r (1 + ^ + (j? ~ ra + mm 
~w=S - v 11 F7i + ,,_ u_ ( - UAL) 

N MM {ki} i=1 i U + Z/ l l ^iJ 

The 2n-point function is given by analytically continuing v = —ifi, changing to MQM 
normalisation using Eq. ( |4.(j| ) (which amounts to removing the power of v in front), 
and finally taking the modulus: 



((T- 1/R T 1/R r) 

with C({ki}) given by Eq. flTXJ). 



ErY/r- TT r (I ~ ^ + ~ ra + Di) 
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